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ABSTRACT 

A new type of quantum simulator is proposed which 
can simulate any quantum many-body system in an 
isomorphic manner. It can actually synthesize a du- 
plicate of the system to be simulated. The isomor- 
phic simulation has the great advantage that the in- 
evitable coupling of the simulator to the environment 
can be fully exploited in simulating thermodynamic 
processes. 

PACS numbers: 89.80.+h, 71.45.-d 



Simulating quantum many-body systems on 
a classical computer is hard in the sense that 
the simulation takes exponentially long time and 
large memory as the size of the system increases 
[1,2]. The difhculty stems from the fact that 
the Hilbert space of the system consists of ex- 
ponentially many states as the function of the 
number of physical variables involved. The ef- 
ficiency of simulating a system by another well 
controlled quantum system (quantum computer) 
was conjectured by Feynman [1] and has been 
justified by Lloyd [2]. Because of the decoher- 
ence problem [3] , despite the recent advancement 
of the quantum computation theory especially on 
quantum error correction [4,5], it remains to see 
whether one can construct a large size quantum 
computer and maintain the delicate quantum co- 
herence in order to accomplish meaningful com- 
putations or simulations, although many efforts 
have been made to implement quantum gates or 
small size quantum computers [6]. 

Here we propose that a universal quantum iso- 
morphic simulator (IS) can be constructed by tai- 
loring the many-body interactions among a suit- 
able physical system. Just for convenience in this 
paper, the system being simulated is called the 
simulatee. One may choose physical systems con- 
sisting of, for example, many isolated quantum 
wires each with a freely moving electron inside, 



or isolated one-dimensional arrays of Josephson 
junctions [7] each having a magnetic fluxon freely 
moving back and forth. The position of the elec- 
tron in the wire or the fluxon in the Josephson 
junction array provides one continuous degree of 
freedom which can be assigned to a continuous 
variable of the simulatee. For simplicity, we will 
use the terms wire and particle in the following, 
whether they are a quantum wire and an elec- 
tron or a Josephson junction array and a fluxon. 
Given the Hamiltonian 

H{t)=H^^H\t), 

where Hq describes the unperturbed system of 
simulatee and H'(t) is the external perturba- 
tion, one can tailor the many-body interactions 
among the system of wires which is called a quan- 
tum isomorphic simulator such that the Hamil- 
tonian describing the simulator is identical (in 
some sense) to H{t), the Hamiltonian of the sim- 
ulatee. This is why it is called an isomorphic 
simulator. First note that Hq — K can be de- 
composed into a sum of two-body interactions in 



general, Hq — K 



Y,i-jf^ijixi^^j)i where Xt,Xj 



are continuous variables which can be simulated 
by two wires (of course, discrete variables can 
also be simulated using wires with double poten- 
tial wells or something else), and the external 
perturbation H' it) is a sum of one-body poten- 
tials, H'{t) = J2k9ki^k), K is the total kinetic 
energy of the simulatee which is naturally mim- 
icked by the inherent inertial motion of the parti- 
cles (electrons or fluxons) of the simulator. One 
may divide each wire into A'^ parts (discretiza- 
tion with precision 1/iV), the position-dependent 
interaction hij{xi,Xj) between two wires can be 
achieved by drawing N'^ connections between the 
two wires so that the interaction potential of the 
ith particle at position Xi and the jth particle 
at position Xj is proportional to hij [xi ,Xj). For 
example, if the wires are Josephson junction ar- 
rays and the particles are magnetic fluxons, one 
may enclose the magnetic fluxes at each pair of 
discrete positions {xi,Xj) by a superconducting 



ring so that the subsystem of the ring and the 
zth and jth fluxons has a total interaction energy 
hij{xi,Xj) when the ith fluxon is at position Xi 
and the jth fluxon is at Xj [8-10]. The one-body 
potentials gk(xk) can be easily realized by ap- 
plying position-dependent external fields to the 
particles. A simulatee with M continuous vari- 
ables can be simulated by an IS consisting of M 
wires, at most ^M{M — 1)N^ connections and 
no more than MN externally applied fields, with 
1/N discretization precision. The cost of phys- 
ical resources is polynomially bounded, and the 
simulator simulates the simulatee in an isomor- 
phic manner, i.e. there is one-to-one correspon- 
dence between the essential Hilbert spaces of the 
two systems, with the eigenenergies of two corre- 
sponding states proportional to each other; the 
two system are governed by the same Hamilto- 
nian (up to a proportion factor) and undergo the 
identical dynamics up to a proportion factor in 
energy and time etc. 

To understand more about what isomorphic 
means, one may compare the IS with the con- 
ventional (quantum) simulator (CS) [2]. On a 
CS, there is generally no one-to-one correspon- 
dence between the eigenstates of the simulator 
and the simulatee with the corresponding en- 
ergies proportional to each other. In particu- 
lar, the ground state of the simulatee is repre- 
sented by the ground state of the IS, while the 
lowest energy state of the CS can not encode 
the ground state of the simulatee because that 
ground state is not known a priori. Although 
a CS may exploit the inevitable coupling with 
the environment to simulate open systems [2], 
the exploitation is strictly limited. For exam- 
ple, when simulating the relaxation process of 
an excited system, a CS should at least preserve 
some energy to encode the information for the 
ground state of the simulatee. These remaining 
excitations make the system unstable and vulner- 
able to perturbations. By contrast, an IS is not 
so niggardly, it happily dissipates all its energy 
and gets to the lowest energy level which rep- 
resents the ground state of the simulatee. This 
example has already demonstrated one scheme 
for an IS to simulate thermodynamic processes. 
Since the IS is just a duplicate of the simulatee 
in the isomorphic sense, its natural relaxation 
naturally mimics the relaxation of the simula- 
tee. This scheme is particularly useful to simu- 



late these relaxation processes where the actual 
forms of the external perturbations are subor- 
dinate, what really important are the internal 
structure of the simulatee and the statistics of 
the perturbations. To this extent, an IS is out- 
standing in simulating the thermodynamics and 
dynamics of complex systems such as biological 
macromolecules [11] and atomic and molecular 
clusters [12]. With a proper Hamiltonian de- 
scribing the interactions among smaller pieces 
of the complex system and some random terms 
mimicking the stochastic perturbations, a clas- 
sical computer may provide some information, 
but a full simulation takes exceedingly long time. 
By contrast, if provided with the correct Hamil- 
tonian, an IS can synthesize a duplicate of the 
complex system, the motion of the duplicate is 
identical to the simulatee in the isomorphic sense 
that there is one-to-one correspondence between 
the eigenstates and the corresponding energies 
are proportional to each other, the dynamics' of 
the two systems are identical up to a constant 
factor in time scale. For instance, if provided 
with the correct Hamiltonian describing the in- 
teractions among the amino acids and their in- 
teractions with water, a synthesized protein on 
an IS should fold itself very fast as the real pro- 
tein does [11]. On an IS, scientists can test their 
models and design new complex systems con- 
veniently. To simulate dynamic processes with 
well-defined driving forces on an IS, dissipation 
should be avoided and one needs to apply exter- 
nal fields with proper forms to simulate the real 
driving forces. 

In conclusion, the proposed quantum simula- 
tor can simulate any quantum many-body sys- 
tem in an isomorphic manner. The cost of phys- 
ical resources is bounded by a polynomial func- 
tion of the size of the system to be simulated and 
the simulation takes a time proportional to the 
time of the real process. The great advantage of 
the isomorphic simulation is that there is a one- 
to-one correspondence between the eigenstates of 
the simulator and the simulatee, the correspond- 
ing energies are proportional to each other, the 
two systems are governed by the same equation 
of motion. In particular, the ground state of the 
simulator corresponds to the ground state of the 
simulatee. Consequently, an IS is outstanding 
in simulating natural thermodynamic phenom- 
ena. The inevitable perturbations from the en- 



vironment can be fully and naturally exploited eds., Nana structures and Mesoscopic Systems (Aca- 

when simulating real thermodynamic processes. demic, San Diego, 1991). 

At present the well-established technology for su- 
perconductive devices is very promising to imple- 
ment an IS using Josephson junction arrays [8- 
10]. The continuing advancement in fabricating 
small size devices may eventually make it possi- 
ble to implement IS' in terms of semiconductor 
nanostructures [13]. It is exciting to expect that 
in the near future quantum isomorphic simula- 
tors are widely used to synthesize and test vari- 
ous complex systems. 
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